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The phase transition of 4D simphcial quantum gravity coupled to U{1) gauge fields is studied using Monte-Carlo 
simulations. The phase transition of the dynamical triangulation model with a vector field {Ny ~ 1) is smooth 
compared with pure gravity (A^v ~ 0). The node susceptibility (x) is studied by the finite size scaling method. 
At the critical point, the node distribution has a sharp peak in contrast to the double peak in pure gravity. From 
the numerical results, we expect that 4D simplicial quantum gravity with 17(1) vector fields has a higher order 
phase transition than 1st order, which means the possibility to take the continuum limit at the critical point. 



1. Introduction 

The phase structure of 4D pure simplicial quan- 
tum gravity has been intensively investigated. In 
4D pure gravity, two distinct phases are known. 
For small values of the bare gravitational cou- 
pling constant the system is in the so-called elon- 
gated phase, which has the characteristics of a 
branched polymer. For large values of the bare 
gravitational coupling constant it is in the so- 
called crumpled phase. Numerically, the phase 
transition between the two phases has been shown 
to be 1st order 0J^. As a result, it is difficult to 
construct a continuum theory. Our next step is 
to investigate an extended model of 4D quantum 
gravity. Recently the phase structure of the ex- 
tended models of 4D quantum gravity has been 
studied numerically [^ji] . In a model with vector 
fields, the intermediate phase has been observed 
between the crumpled phase and the elongated 
phase. We consider the possibility of a continuum 
limit at the critical point. In order to investigate 
the phase transition, we measure the finite size 
scaling of node susceptibility (x) , 

X = (< A^o' > - < iV o >')M, (1) 
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as well as the scaling property of the mother 
boundary in the case of 4D simplicial quantum 
gravity coupled to one gauge field {Ny — 1). 

2. Phase Diagram with Gauge Fields 

We consider the partition function of simplicial 
gravity coupled to U{1) gauge fields. The total 
action is S* = Seh + Spi. We use the Einstein- 
Hilbert action for gravity, 

Seh[A, G] = Jd^x^iA - ^R), (2) 

where A is the cosmological constant and G is 
Newton's constant. We use discretized action for 
gravity, 

Seh[k2,K4] = KiNi - K2N2, (3) 

where K2 ^ K4 is related to A and Ni is 

the number of i-simplices. We use the plaquette 
action for U{1) gauge fields ||^, 

Spi = ^o(t.,fe)[A(/,,) + A{l,k) + Aih^)]^ (4) 

where hj denotes a link between vertices i and 
j, tijk denotes a triangle with vertices i, j and 



k, A{lij) denotes the U{1) gauge field on a link 
lij, and o{tijk) denotes the number of 4-simplices 
sharing triangle i^fc- We consider that a partition 
function of gravity with Ny copies of the U{1) 
gauge fields is 

^4 t(2_D)eT(4_D) 

n/ n dA{l)e-''^K (5) 

Nv l£t{2D) 

We sum over all 4D simplicial triangulation, 
T(4:D), in order to carry out a path integral over 
the metric. Here, we fix the topology with S'^. 
Numerically, in the case of adding vector fields, 
three phases have been found |^,||. A schematic 
phase diagram has been shown in ref.|Q. An in- 
termediate region is called the smooth phase be- 
tween these two transition points,]] and We 
expect that the phase transition at Kg is continu- 
ous and leads to continuum limit of 4D quantum 
gravity. Now, let us notice the transition at Kg in 
the case of Ny = 1- 

3. Numerical Analysis of Phase Transition 

(Nv = 1) 

In this section we report on two numerical ob- 
servations: the node susceptibility (x) and the 
histogram of A*o. In Fig.l we plot the node sus- 
ceptibility (x) as a function of K2 with volume 
N4 = 16K, 2AK and ?,2K, respectively. The node 
susceptibility (x) has a peak value at the crit- 
ical point (Kg)- We measure the peak value in 
each size. As a finite size scaling, the peak value 
[Xmax) and the width of peak (^Kg) grow as iV4 
in power. The susceptibility exponents, A and F, 
are defined by B 



Xmax OC ((5k2 OC A'"4 ^). 



(6) 



From the numerical result (Fig.l), we obtain the 
susceptibility exponent A = 0.4(1), (F - 0.5(3)). 
These values are apparently smaller than 1. In 
Fig. 2 we show the histogram of Nq at the crit- 
ical point (Kg — 1.37147(1)) with the size of 



X N4=16K 
+ N4=24K 
% N.=32K 
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Figure 1. Node susceptibility (x) plotted ver- 
sus the coupling constant (K2) and Xmax plotted 
versus A^4 with log-log scales for Ny = 1. 

A'4 — i2K. In pure gravity, a double peak struc- 
ture has been found||^J^. The fact shows that the 
phase transition is 1st order. However, in the case 
of Nv — 1, the double peak structure disappears. 
We then consider that the phase transition be- 
tween the crumpled phase and the smooth phase 
may be continuous, not 1st order. 

4. Scaling Property of 4D Simplicial Quan- 
tum Gravity 

In this section we discuss the scaling structure 
of 4D DT manifolds, focusing on the scaling struc- 
ture of boundaries in 4D Euclidean space-time us- 
ing the concept of geodesic distances. In order to 
discuss the universality of the scaling relations, 
we assume that the boundary volume distribu- 
tion {p{x, D)) is a function of a scaling variable, 
X — V/D"-, with the scaling parameter a in the 
analogy of 2D quantum gravity. Here, V denotes 
the volume of the boundary and D is the geodesic 
distance. The expectation value of the boundary 
three-dimensional volume appearing at distance 
D has been introduced in ref.lHl, 



< 



1 r dVVpix = V/D'^,D), 

Jvn 



(7) 



^The usual phase transition point (kj) and jj^e obscure 
phase transition point (kj) a^r^ defined in ref. [ 



where vq denotes UV cut-off of the boundary vol- 
ume and N is the normalization factor. If the 
boundary volume has a scaling property with the 
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Histogram of nodc(A^o) at the critical 
= 1.37147) with the size of = 32K. 



universal distribution p(x, D) and wq 0, 



(8) 



Then, we obtain a finite fractal dimension, rf/ = 
a + 1, with the fractal dimension df. We mea- 
sure the volume of the mother boundary as a 
function of D. The mother boundary is defined 
by the boundary having the largest tip volume. 
In Fig. 3 we plot the mother boundary volume, 
< V^'^^ >, at the critical point with a size of 
Ni — 32K. As a result, the mother boundary 
volume shows scaling and we obtain the scaling 
parameter a — 3.7(5). Then, we can estimate 
the fractal dimension [df = 4.7(5)). On the 
other hand, we measure the Hausdorff dimension, 
dn = 4.6(2). Both results are consistent. Thus, 
we expect that the boundary volume distribution 
has a scaling property in the sense of the manifold 
at different distances from a given 4-simplex and 
looks exactly the same after a proper rescaling of 
the boundary volume. 

5. Summary and Discussions 

Let us summarize the main points made in the 
previous sections. For the phase transition, we 
show the finite size scaling at the critical point 
and the histogram of node. From the numeri- 
cal results, the phase transition is smooth in con- 
trast to pure gravity. Also we show the scaling 
property of the mother boundary, where the scal- 
ing parameter is consistent with the fractal di- 
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Figure 3. Power fit to mother boundary three- 
dimensional volume < T^^^^ > at the critical point 
(k2 = 1.37147) with the size of iV4 = 32i^ with 
log-log scales. 

mension. We expect that the boundaries have a 
fractal structure and the universality of the scal- 
ing relations. From a modification of the Balls- 
in-Boxes model we expect that the simplicial 
quantum gravity coupled to matter fields will 
have the possibility of a continuous phase transi- 
tion in Nv > 1. We expect that the existence of 
genuine 4D quantum gravity at the critical point, 
^2, with the vector fields. 
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